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Abstract 

Universal quantum cloning machines (UQCMs), sometimes called 
' quantum doners, generate many outputs with identical density matri- 

, ces, with as close a resemblance to the input state as is allowed by the 

t-H ' basic principles of quantum mechanics. Any experimental realization 

, of a quantum doner has to cope with the effects of decoherence which 

terminate the coherent evolution demanded by a UQCM. We examine 
how many clones can be generated within a decoherence time. We 
compare the time that a quantum doner implemented with trapped 
ions requires to produce M copies from N identical pure state inputs 
and the decoherence time during which the probability of spontaneous 
emission becomes non-negligible. We find a method to construct an 
^ N — » M cloning circuit, and estimate the number of elementary logic 

£ , gates required. It turns out that our circuit is highly vulnerable to 

spontaneous emission as the number of gates in the circuit is exponen- 
, tial with respect to the number of qubits involved. 

•i-H . 

1 Introduction 

Although the "no-cloning theorem" by Wootters and Zurek Q states that 
it is impossible to copy arbitrary quantum information perfectly, quantum 
mechanics does allow us to make approximate copies. Since Buzek and 
Hillery first presented a unitary transformation, known as a universal quan- 
tum cloning machine (UQCM), to make two identical approximate copies 
of an input qubit Q, the subject of quantum cloning has been investigated 
intensely (for example, see [||, ||, |(|) and experiments have recently been 
carried out [Q, ^, [| . 

One application of quantum cloning is in quantum computation, as copy- 
ing is a fundamental process in information processing. It has been shown 
that there are some quantum-computational tasks whose performance can 
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be enhanced by making use of quantum cloning |1C]. However, physical 



realizations of quantum circuits are always fraught with difficulties due to 
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decoherence, as it is virtually impossible to isolate a quantum system from 
its environment perfectly. 

Here, we investigate the circuit complexity of universal quantum cloning 
and estimate the time T a quantum cloning machine requires to make M 
copies of N identical pure state inputs. By comparing T with the decoher- 
ence time 7dec of the qubits in a physical system, we can estimate how M 
or N will be restricted and how much quantum information can be copied 
practically. In order to estimate the decoherence time, one needs to model a 
specific realization, and several are possible: cavity QED, parametric down- 
conversion, and ion traps. To provide a concrete example, we examine an 
ion-trap realization of the UQCM network, and assume all instrumental 



sources of decoherence [11, |12|, |13|] have been eliminated, leaving us with 
spontaneous emission as a well-characterizable source of decoherence. 

The analysis of the times T and Td ec is analogous to that reported in 
[14 1, in which upper bounds are determined for the bit size L of the number 
to be factorized by using Shor's algorithm. 

In the next section, we briefly review the ideas of quantum cloning, 
and then, in Section |3[ count the number of elementary logical operations 
needed to build a quantum circuit for N — > M cloning. In Section ||, we 
compare the cloning time T with the decoherence time Td ec of the quantum 
circuit when it is implemented with trapped ions, as in [14]. Other sources 
of decoherence (vibrational quanta damping, heating, inhomogeneous trap 
fields, and laser fluctuations) will also lead to limitations in the number of 
good quality copies such a UQCM can generate. 



2 Quantum Cloning Transformation 

The task of a general universal quantum cloning machine is to copy N 
identical pure states, described as p m = 

(\ip)(ip\) m , into M-particle output 
states (M > N), p ont , with the following conditions: 

i. ) The reduced density operators for any one of the M outputs are iden- 

tical to each other, i.e., 

P° ut = Pi Ut > (1) 
where p° ut is a reduced density operator with respect to the i-th par- 
ticle. 

ii. ) The quality of the copies does not depend on the input states, i.e. the 

fidelity between the input and the output states is independent of the 
input state. The fidelity F is defined by F = (ip\p° nt \^) . The word 
"universal" refers to this condition. 

iii. ) The copies are as close as possible to the input state as a natural 

requirement for a cloning machine. Thus, the fidelity F should be as 
close as possible to 1. 
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Here, we consider only a two-dimensional system (qubit) as a physical 
system, as systems of higher dimension are beyond the scope of current 
candidates for the realizations of quantum circuits. 

The cloning transformation of Buzek and Hillery ||] makes two copies 
from one original qubit. It is written as 

\o) a \o) b \o) x - y||o) |o) 6 |o) ie + y|(|o) |i) 6 + |i) |o> 6 )H> ai 

\l)a\0) b \0)x -> \/||l>o|l>6|l>* + y|(|0)«|l)6+|l>«|0> 6 )|0) a! . (2) 

In Eq. (|2j), the first qubit with subscript a is the state to be copied, the 
second one labelled b is the blank paper that becomes one of the copies after 
the process, and the qubit with x is an ancilla bit which can be regarded as 
the state of the machine. The fidelity of this process is found to be |, which 
is independent of the input state, as desired. 

Buzek et al. also presented a way to construct a quantum network for 
this UQCM ml (Figure [j](a)). In Figure [j](a), R is a single qubit gate which 
rotates the basis vectors by an angle 9 as 

R(8)\0) = cos#|0) +sin0|l>, 

R(0)\1) = -sin0|O> + cos0|l>, (3) 

and • and © symbols connected with a line denote a controlled NOT gate 
(CNOT) with • and © as a control bit and a target bit, respectively. By 
adjusting the rotation angles of three single qubit gates, any two of the three 
states at the output represent copies of the qubit a. 

Eq. (||) was generalized by Gisin and Massar |J to produce M copies 
out of N inputs. This transformation is described by 

M—N 

u N , M m)\o)® 2(M ~ N) = Y. a ^ M -^^^)®Rm^ (4) 

3=0 

/N + l l {M -N)\(M - j)\ 
ttj ~ V MTTy (M-N-j)\M\ ' 

RjW = \(M — N - j)ip*,j(ip*) ± ), 

where |iV^) is the input state consisting of N qubits, all in the state \ip). 
\{M — j)if), jip^) is the symmetric and normalized state with M — j qubits 
in the state \ip) and j qubits in the orthogonal state |V'" L )- Rji^) represent 
the internal state of the machine and Rj(ip) -L Rk{ip) holds for all j 7^ k. 

3 Generic Quantum Cloning Circuit 

Buzek et al. constructed a quantum circuit for 1 — ► M UQCM which explic- 
itly realizes Eq. (0) by generalizing the corresponding circuit for the 1 — ► 2 
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Figure 1: Quantum circuits for (a) 1 -> 2 UQCM, and (b) 1 -> M UQCM. 
These circuits consist of two stages, namely the preparation and the copying 
stages. The preparation stage provides appropriate amplitudes for orthonor- 
mal bases according to Eqs. (2) or (4). The cloning stage permutes those 
amplitudes among all qubits, entangling them([13, 14]). 
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Figure 2: Quantum circuit for N — > M UQCM? It turns out that this cloning 
circuit does work properly for any M when N = 1, N = 2, and many other 
N and M as long as the condition (|6|) is fulfilled. If we are allowed to 
introduce additional auxiliary qubits inside the doner, this circuit works for 
all combinations of N and M(> N). 



case [16]. Figure |l](b) shows their circuit for 1 — > M cloning machine. It is 
basically a natural extension of the 1 — * 2 circuit, consisting a preparation 
stage and a cloning (permutation) stage. Let us consider if it is possible to 
extend this 1 — > M cloning circuit further to an N — > M cloning circuit, 
by making up a circuit in Figure ^ in which the preparation stage is built 
up with both single and two qubit operations and the cloning stage is a 
sequence of CNOT or similar multi-qubit operation gates. 

Eq. (H) gives us a few hints on the possible construction of the quantum 
circuit for an N —* M UQCM. First, M — 1 qubits are necessary to represent 
the internal states of the machine. The total number of qubits needed to 
implement this transformation is 2M — N. 

Second, every basis of the form of \(M — jV'~ L ) or Rjiip) i n Eq. @ is 
a basis in the symmetrical subspace of the Hilbert space 
where 7i is the space spanned by |0) and |1). Although the symmetrical 
subspace is a rather small subspace in the whole Hilbert space, almost all 
computational bases are involved in Eq. (|]), especially when iV is small. 

Therefore, almost 2 2M ~ N amplitudes need to be distributed to appro- 
priate computational bases to perform the transformation (B), although the 
number of distinct amplitudes, ay, is only M — N + 1. The fact that such a 
great number of amplitudes are necessary imposes a certain condition upon 
the values of N and M in order for an iV — > M cloning machine in Figure § 
to work properly. 
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Because the task of the cloning stage which consists of a sequence of 
CNOTs is to permute the amplitudes of all basis vectors which span the 
whole Hilbert space, all the amplitudes in Eq. (||), aj, must be provided by 
the preparation stage. Therefore, the preparation stage should be able to 
generate an entangled 2(M— iV)-qubit state with a set of any real amplitudes 
that may appear in Eq. (Q) by adjusting the rotation angles of each single 
qubit gate. 

For example, in the 1 — ► 2 UQCM (Eq. (|2|) and Figure ffl(a)), the am- 
plitudes appearing in the output states are y|, y^t, and Thus, the 
preparation stage needs to provide these three coefficients and it turns out 
that, with the configuration of the cloning stage of Figure [j](a), the two-qubit 
state emerging from the preparation stage should be 

|^)p-p = y|| 00 ) + yi| l) + y^lll), (5) 

when the third qubit is used as the internal state of the machine. The rota- 
tion angles of the single qubit gates in this stage are determined accordingly. 

However, this scheme only works when the number of distinct computa- 
tional bases in Eq. (||) is smaller than the dimension of the Hilbert space 
that the preparation stage deals with. Otherwise, the preparation stage 
cannot generate enough amplitudes required by Eq. (||). This condition is 
written as 

where the LHS represents the number of bases that appear in Eq. (0) and the 
RHS is the dimension of the Hilbert space where qubits in the preparation 
stage lie. 

Eq. (JgJ) is fulfilled for any M when N = 1, which is the case of Figure 
|(b), and when N = 2. Still, it becomes rather complicated when it comes 
to other combinations of N and M. An example in which we can see the 
violation of Eq. (||) clearly is the case of 2N = M, where the LHS of Eq. 
(£]) for N S> 1 can be estimated as 




(7) 



This always exceeds the corresponding RHS of Eq. (||), 2 2N . Nevertheless, 
Eq. (^) is satisfied for M which are sufficiently large compared with N, 
since its LHS is approximately 2 2M ~ l /y/irM , which is smaller than the RHS 
when M > vr"^ 4 ^- 2 . 

Note that the condition (^) is not necessary if we are allowed to make 
use of more auxiliary qubits inside a cloning machine. Providing 2M — N 
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qubits at most, instead of 2(M — N), enables the preparation stage to gen- 
erate enough number of amplitudes for the cloning transformation. Then 
the cloning stage can complete the whole process by allocating those am- 
plitudes to appropriate bases, leaving the auxiliary qubits disentangled, in 
state |0), from the legitimate 2M — N output qubits. Introducing these aux- 
iliary qubits does not affect the estimation of the number of gates, which is 
discussed in the following. 

Let us now count the number of elementary gates in the quantum cloning 
circuit, especially CNOT gates, as a CNOT gate takes a longer time to be 
performed than a single qubit gate 14]. 

The preparation stage should be designed so that it generates an arbi- 
trary set of real amplitudes whose number is written by the LHS of Eq. (||) . 
As far as the quantum cloning transformation Eq. (||) is concerned, the am- 
plitudes should not necessarily be arbitrary as they are specifically described 
as in Eq. @). However, in order to keep the full "controllability" on the 
output states, as in the case of 1 — > 2 UQCM, we assume that the prepa- 
ration stage can generate arbitrary superpositions of 2 dprcp bases with real 
amplitudes, where d prep = 2 2 ( M ~ N ^ is the dimension of the 2(M — A)-qubit 
Hilbert space. 

The transformation that the preparation stage performs can be written 
as 1 00 • • • 0) — ► Yl'kZ'o 1 c k\k), where are real numbers complying with the 
normalization condition ^2 k c| = 1. This is a unitary transformation which 
can be described by a d prep x d pTep matrix, on the d prep -dimensional initial 
state vector (1, 0, • • • , 0) T , and its components are given by u m i = c m and 
the rest of them are arbitrary for our purpose. 

It is known that an outright circuit implementation of a d x d unitary 
matrix requires, in general, 0(d 2 (log d) 2 ) elementary operations 17]. How- 
ever, a more efficient circuit to create an arbitrary quantum superposition 
starting from 1 00 • • • 0) has been proposed in |l8| and its complexity is given 



by 0(d(log d) 2 ). Thus we simply take 0(d prep (log d prep ) 2 ) as the number of 
CNOTs in the preparation stage of a UQCM in the following calculations. 

The cloning stage is also generically hard to construct. The only excep- 
tion we know of is the case of 1 — > M UQCM, which can be realized by 
the circuit in Figure ||(b). In a more generic N — > M case, we can build 
up a circuit as follows. As mentioned above, the cloning stage only per- 
mutes all bases with non-zero amplitudes among all computational bases 
of n®{ 2M - N ) _ Thus, we can estimate the number of gates with only the 
knowledge of number of bases involved, even if we have no information on 
the actual permutation, i.e. which basis goes to which. 

Let us take the 1 — > 2 quantum cloning circuit as an example to simplify 
our description of its construction, although its efficient circuit is already 
given in Figure |l[ Table |l| shows all the necessary permutations of bases 
to complete the transformation from Eq. (|5]) to Eq. (0). This permutation 
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Table 1: The permutation of computational bases performed in the cloning 
stage of the 1 — > 2 quantum cloning circuit. The choice of bases on the left 
hand side follows Eq. (5). 
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Figure 3: The quantum circuit for the permutation of bases. This figure 
depicts a circuit for the cloning stage of 1 — > 2 UQCM. An additional ancilla 
qubit is introduced as a flag to flip bits. Although there exists a much more 
efficient circuit for the same process, which is shown in Figure 0(a), this 
circuit can be applied to the UQCM with more input qubits, as in Figure ||. 



can be implemented as a circuit by carrying out each transformation one by 
one. Since 1 010) is not used in the preparation state, 1 101) — > 1 010) should 
be performed first, otherwise any other destination state may be used as 
a preparation state later and thus the preceding transformation comes to 
naught. A basis which is exempted from being used in the preparation, 
such as 1 010) in this case, always exists in a more general iV — > M case 
as the equality in condition (^) never holds for N, M E IN. Hence, the 
appropriate order of transformation is 1 101) — * 1 010) , |001) — * 1 101) , 1 1 1 1) — * 
|001),|100) -> |111),|001) -> |100), where |001) is used as a buffer in the 
last three operations, as it is not used in the final state, to swap 1 100) and 
|111). 

Figure |3| shows a general way to permute bases. In this figure, an unfilled 
circle, o, denotes a control bit in-between two NOT gates. Unlike a filled 
circle •, o activates connected operations when the bit value at o is 0. 

Here, an ancilla qubit is introduced in order to flip the incoming binary 
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numbers according to the permutation required for a UQCM, as shown in 
Table |l|. Finding a circuit for the basis permutation without ancillas is gen- 
erally very hard, hence we make use of it here for convenience. The ancilla 
is prepared by the UQCM in a state |0) and it ends up in the same state 
|0) after the whole process. As the ancilla qubit is necessarily disentangled 
by the end of the process, we do not have to take it as a part of the output 
state. 

In the earlier paragraph of this section, we have identified the generic 
features required of a copying circuit. For concreteness, it is important to 
show a specific example and in the following we do this for a 2 — > 4 doner. 
Our circuit is not unique and other examples of course can be found. 

The transformation of the state 1 00) by the 2 — > 4 quantum cloning can 
be written from Eq. (0) as 

|oo)|oooo) -> y?|oooo)|ii) 

+ J Jr(|0001) + |0010) + |0100) + |1000»(|01) + |10», 
V oO 

6 bases 



+y-(|0011) + --- + |1100))|00), (8) 

and the transformation for 1 11) 1 0000) can be obtained by flipping every bit 
in Eq. (||). 

As mentioned earlier in this section, we need four qubits for the prepa- 
ration stage in addition to two input qubits, thus, six qubits in total. The 
preparation stage generates all amplitudes for the 15 bases appearing in Eq. 
(||). The output from the preparation stage is, for example, 

8 bases 



|0000) -► |^) prep = W -10000) + \ — (|0001) + ••• + |1000)) 
V 5 V oO 



6 bases 



+y_(|1001) + --- + |1110)). (9) 

The quantum circuit to carry out this transformation is given in [[03] and is 
also depicted in Figure f|. The boxes represent a single qubit operation Ug 
of the form 

V sin a — cos a J 
and all 6 in Figure |I| are given by 

ho = tan" 1 J^, 02i = tan" 1 J^, 
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Figure 4: The preparation stage for 2 — > 4 quantum cloning circuit. Angles 



6 in boxes are given by Eqs. (11). 
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One possible permutation of bases that achieves the transformation (||) 
from |-i/;) pre P in Eq. (||) is shown in Table ^. This can be implemented by the 
quantum circuit in Figure [|. With the quantum circuits of Figures |] and 
[|, the quantum cloning transformation (||) is carried out faithfully and the 
upper three qubits in Figure || will be the clones, while the next two qubits 
are the state of the doner and the lowest disentangled qubit is an ancilla 
which can be discarded. 

We can now estimate the number of CNOTs in the whole cloning cir- 
cuit. In p9| , the number of basic operations, i.e. single qubit operations 
and CNOTs, to simulate a multi-qubit controlled operation has been given 
as 0(n 2 ) for a controlled-?/ gate with n — 1 control qubits and one auxil- 
iary qubit. Therefore, the upper bound for the number of CNOTs in the 
cloning stage is the order of the number of bases involved, which is ap- 
proximately twice of the number of amplitudes produced in the preparation 
stage, multiplied by the square of the number of qubits. As the number 
of the bases is given by the RHS of Eq. @, 0(2 2M M~2), and there are 
0(2 2 ( M ~ N ' (M — N) 2 ) CNOTs in the preparation stage, a quantum cloning 
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Table 2: The permutations of bases to complete the transformation Eq. (||) 
when the circuit of the preparation is given by Figure |j. Many of them are 
omitted in this table since they are almost obvious by comparing Eq. (||) 
and Eq. @. 



Initial basis 


-> Final basis 


Initial basis 


-* Final basis 


1 000000) 


-> |000011) 


1 110000) 


-> |111100) 


000001) 


-> 000101) 


110001) 


111001) 


000010) 


-> 000110) 


110010) 


111010) 


000011) 


-> 001001) 


110011) 


-> 110101) 


|001101) 


-> 1 101000) 


|111101) 


-» |101011) 


001110) 


-» 110000) 


111110) 


-> 110011) 




Figure 5: The cloning stage for the 2 — * 4 quantum cloning circuit. Only 
four permutations |000011) |001001), |000100) |001010), |000000) 
|000011), and 1 110011) — > 1 110101) are shown. Each of these permutations 
consists of three multi-qubit control operations. This figure represents only 
one possible example because the order of permutation is not unique. 



11 



circuit of the type of Figure | contains 0{2 2M (M - N) 2 (2- 2N + M~a)) 
CNOTs at most in total. 

This circuit looks rather inefficient especially because of the preparation 
stage. The inefficiency comes partly from our requirement that the prepa- 
ration stage should be able to generate arbitrary superpositions. If both N 
and M can be fixed and we do not have to control the parameters of each 
single qubit gate, the task of the preparation stage is much easier and the 
configuration of its circuit may well be much simpler and more efficient. 

The cloning stage is also inefficient due to the complexity of the network 
for a general permutation. Despite the linear dependence on M in the case of 
1 — ► M UQCM, the number of CNOTs grows exponentially as M increases 
when there are iV inputs, since each permutation is done one by one in our 
circuit. It might be possible to find a more efficient circuit, however, we 
have not succeeded. 



4 Decoherence Time and Cloning Time 

With the number of CNOT gates estimated in the previous section, we can 
now compare the cloning time T with the decoherence time Td ec . We focus 
on a quantum computing realization which makes use of cooled trapped ions 



[ 20 1 . In the following, we assume that spontaneous emissions are the only 
source of decoherence, and we only discuss the process without error correc- 
tion codes. Going beyond these constraints will be discussed elsewhere. 

The Hamiltonian operator for a two-level ion of mass m, interacting with 
a phonon as a result of the centre-of-mass (cm.) motion with frequency u, 

" = ^]#=wt[i i>(0| ' 1 +i o)(1 I'' , i> (12) 



is then given by [2C, 14], 



where rj = (2n/ A) ^Jh/2mv is the Lamb-Dicke parameter, Cli is the Rabi 
frequency of the <-» 1 transition with and 1 denoting the ground and the 
excited states of the ion. The a and a) are the annihilation and creation 
operators of the phonon. As in the original proposal |2C[] , we assume that 
qubits are encoded in the internal state of ions and the phonons are used as 
the information bus. The denominator y/2M — N is a consequence of the 
fact that an N -> M UQCM network needs 2M - N qubits (ions). 

The elementary time step for a CNOT gate with this system can be 
written as 

(13) 

The total processing time for cloning is 
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Figure 6: Three-level model of the ions used in a quantum cloner. The 
transition between and 1 represents the qubit. An external laser drives 
this transition with Rabi frequency fix. The laser inevitably couples level 
to other non-resonant levels such as level 2. The Rabi frequency of this 
transition is and level 2 decays with a rate IV 



where e is a some proportion factor. Here, we use the number of CNOTs in 
the preparation stage, because it is usually dominant over that of CNOTs in 
the cloning stage, as mentioned above. Also, we assume that all the CNOTs 
are performed sequentially one by one, despite the possibility of building up 
a circuit which performs several operations in parallel. Thus, the discussion 
below gives only a naive upper bound for the number of cloneable qubits. 
To minimize T, we wish to increase the value for which is related to the 
decay rate of the excited state by [14, ^T] 



67rc 3 e 2 



E\ 



(15) 



where E,c,eo, and u>i are the electric field strength of the laser, the speed 
of light, the permittivity of vacuum, and the transition frequency between 
the states and 1, respectively. Clearly, we would like to minimize the 
decoherence effects of Ti by going to a metastable level, but then Qi is also 
small. We can increase the Rabi frequency by increasing the electric field 
strength of the laser E, but it would then cause transitions to other higher 
levels which may be more rapidly decaying than state 1, and eventually 
spontaneous emissions which destroy the coherence of the system. Let |2) 
represent all other auxiliary levels which are coupled to the ground state |0). 

As a stronger laser increases the population in the auxiliary level, it 
increases the rate of spontaneous emissions from level 2. Thus, we compute 
the probability of an emission from either level 1 and 2, and then minimize 
this probability to have an intensity-independent limit to the number of 



output qubits under the effect of spontaneous emission [21]. 

The probability of a spontaneous emission from the upper level of the 
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qubit during the cloning process is 



Pl ^ = -2ri(2M-Jv)r. 



(16) 



The factor | is present because we can assume that on average half of the 
qubits are in the upper level during the whole cloning period T. Because 
the auxiliary level is populated only when interacting with the laser, the 
probability of a spontaneous emission from level 2 is written as 



n2 



(17) 



where A2 is the detuning between the frequency of the laser and that of 
the transition <-> 1. We now obtain the total probability of spontaneous 
emission 



Ptotal = Pl^O + P2^0 

47re v / (2M - N)*T lt 
V 
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where x = fii/\/Ti and we used 
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2 2N 



(18) 



W2V n| 



which is derived from Eq. (J15|) . The minimum value of ptotal with respect 
to x is 



P mi n = - ( ^ ) 5 ^2^(2M - iV )( M - iV) 2 ( J_ + 1 



(19) 



In order for the cloning process to be successful, p m \ n should be very small 
compared with unity. This requirement gives the upper bound on M for a 
given N in the expression of 



: (2M n 



N)(M n 



iV) 2 



1 



2 2N 



+ 



1 
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3 



2 4 ^e V^J r 2 



(20) 



Table y shows the results of the actual calculations of the RHS of Eq. (|20|) 
for those ions often utilized in trap experiments along with their numerical 
data. Assuming e= 100 (2^,77 = 0.01 and A 2 = 10 13 [s" 1 }, we find that Eq. 
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Table 3: Atomic data of several possible systems for a quantum cloning 
machine and the RHS of Eq. (|20|), which should be smaller than the LHS. 



We can see that Eq. (20) has no solutions for N > 1 and M > 2, as its LHS, 
which is equal to 31.15 when N = 1 and M = 2, is always larger than the 
RHS shown in this table. Atomic data are taken from [21]. 



Ion 


Ca+ 


Hg+ 


Ba+ 


level 
level 1 
level 2 


4*^1/2 

3d 2 D 5/2 


5d w 6s 2 2 S 1/2 
5d 9 6s 2 2 D 5/2 

5d l0 6p 2 2p i/2 


Gs 2 S 1/2 
5d 2 D 5/2 
Qs 2 P 3/2 


Is''] 

L02 IS' 1 ] 

r 2 [s- 1 ] 


2.62 x 10 15 
4.76 x 10 15 
67.5 x 10 6 


6.7 x 10 i! > 
11.4 x 10 15 
5.26 x 10 8 


1.07 x 10 i! > 
4.14 x 10 15 
58.8 x 10 6 


RHS of Eq. @ (rj = 0.01) 


0.72 


0.084 


2.58 



(20) has no solutions for M > 2 as the LHS of Eq. (f20|) is equal to 31.2 when 
N = 1 and M = 2. However, if we consider only the 1 — > 2 cloning network, 
where the number of CNOTs is six according to Figure [l]( a) , the probability 
of spontaneous emission during the whole process becomes rather small. 
For example, p m in = 

0.062 for Ca+ and 

Pmin — 0.017 for Ba~*~, therefore the 

cloning may be successful. 

Overall, we see that even for a small number of outputs copied from one 
input the decoherence due to spontaneous emissions plays a critical role. If 
we make an optimistic assumption for rj, as rj = 1.0, 1 — ► 2 cloning with 
Ca + , 1 — > 2 and 2^3 with Ba + would become possible. 

If we are allowed to have one more auxiliary qubit, the number of basic 
operations to simulate a (n — 2)-controlled operation can be reduced from 
0(n 2 ) to 0{n) [19]. Replacing the corresponding factors in the above cal- 
culations shows that, with r/ = 1.0, 1— >2, 1— >3, and 2 — ► 3 cloning may 
be possible with Ba + as well as 1 — ► 2 and 1 — > 3 cloning with Ca + . 



5 Summary 

We have first investigated a possible method to construct an N — > M quan- 
tum cloning circuit in order to estimate the number of CNOT gates in it. We 
have estimated the number of CNOTs as 0(2 2M (M - N) 2 {2~ 2N + M~^)), 
with N and M the numbers of input and output qubits, respectively. There- 
fore, the number of gates in the quantum cloning circuit of the type of Figure 
||(b) or Figure ||] is always exponential with respect to the number of output 
qubits. 

With the circuit complexity we obtained, it has been shown that the 
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quantum cloning may be implemented by using the system of trapped ions 
for only a few combinations of small N and M, when the system is not im- 
mune to decoherence due to spontaneous emissions, provided sophisticated 
quantum error correction codes are not used. As spontaneous emissions 
are the only source of decoherence we have considered here, the number 
of possible combinations of N and M might be lowered further by taking 
into account of other effects, such as phonon decoherence |23|], the random 
phase fluctuations of the lasers, the heating of the ions' vibrational motion 
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Nevertheless, unlike the case of factorization by Shor's algorithm 14], 
producing many clones is not necessarily what we expect from the quantum 
cloning. Since even a few copied qubits may be useful in quantum informa- 



tion processing [10], our results should not be interpreted too pessimistically. 
Furthermore, the use of quantum error correction codes will surely ease the 
condition for upper bounds. 



In [25], it was shown that a quantum information distributor, which is 
a modification of the quantum cloning circuit, can be used as a universal 
programmable quantum processor in a probabilistic regime. From this point 
of view, our estimation on the upper bound for the clones implies physical 
bounds on the realization of such a processor with trapped ions. 

One interesting subject we have not considered here is the effect of deco- 
herence on the quality of the clones. If we allow a processing time which is 
longer than the decoherence time, the fidelity between the input and output 
states will surely be lower than what we expect from Eq. (Q). We might 
be able to find a tradeoff between the fidelity and the number of cloneable 
qubits to have "useful" clones in terms of practical quantum information 
processing. 
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